A complete criterion for separability detection 



OS 

o 
o 

(N 

in 



^ ■ 

a,; 

^— > ■ 

c ■ 

c3 ■ 
=5 

cr 



> 
in 
m 

(N 

o 
o> 
o 



X 



Miguel Navascues 1 , Masaki Owari 1,2 and Martin B. Plenio 1 ' 2 
1 Institute for Mathematical Sciences, Imperial College London, SW7 2PE, UK, 
2 QOLS, The Blackett Laboratory, Imperial College London, Prince Consort Rd., SW7 2BW, UK 

Using new results on the separability properties of bosonic systems, we provide a new complete 
criterion for separability. This criterion aims at characterizing the set of separable states from the 
inside by means of a sequence of efficiently soivable semidefinite programs. We appiy this method 
to derive arbitrarily good approximations to the optimal measure-and-prepare strategy in generic 
state estimation problems. Finally, we report its performance in combination with the criterion 
developed by Doherty et al. [l|] for the calculation of the entanglement robustness of a relevant 
family of quantum states whose separability properties were unknown. 

PACS numbers: 



Research on separability criteria, that is, on compu- 
tational methods to determine whether a given state is 
separable or entangled, is a popular subject in Quan- 
tum Information Theory. Starting from the famous PPT 
(Positive Partial Transpose) criterion, a considerable 
number of different separability criteria have been dis- 
covered (see the references in [1, 0; S Q ) ■ Unfortunately, 
the most efficient ones happen to be partial criteria, in 
the sense that they only detect entanglement or separa- 
bility in certain situations. This is not surprising, in view 
of the fact that the separability problem is NP-hard [7| , 
1- 

A complete criterion for entanglement detection is an 
algorithm or method that allows to characterize the set of 
separable states with arbitrary precision. Alternatively, 
we may say that a complete criterion for separability can 
solve any instance of the Weak Membership Problem of 
separability ((J, i.e., the problem of determining if a given 
quantum state is close to the core of the set of separable 
states. Up to now, there exist several different complete 
separability criteria [2, 0,0, 11, 13, 13, 14 1. In all these 
methods, the set of separable states is approximated suc- 
cessively by an appropriate sequence of sets of states. 
The complexity of characterizing each of these sets in- 
creases as we move along the sequence, that converges to 
the set of separable states in the asymptotic limit. That 
way, we can characterize the set of separable states up 
to a precision that is only limited by our computational 
resources. 

While traditionally Entanglement Theory has focused 
most of its efforts in developing criteria for entanglement 
detection [l|, [l2|, EH , thus characterizing the set of sep- 
arable states from the outside, there exist a few complete 
criteria that try to approximate the set of separable states 
from the inside, i.e., whose aim is to detect separability 
instead of entanglement. Examples of the latter type are 
the algorithms invented by Hulpke et al. [l(|, Brandao 
ct al. [11] or Spedaiieri 13]. The main drawbacks of the 
first two are their high time complexity, as estimated in 
The method designed by Spedaiieri, although very 
promising, does not currently have any proven bounds 



on its speed of convergence, and cannot be extended to 
deal with multipartite entanglement [l3j |. 

Let us also remark that, in all the above cases, the main 
algorithm works by almost directly invoking the defini- 
tion of separability, i.e., no insight from Entanglement 
Theory itself was employed in their conception. 

In this paper, we present a new complete criterion for 
separability which takes inspiration from the symmet- 
ric extension criterion developed by Doherty, Parrilo and 
Spedaiieri (the DPS criterion) However, while the 

DPS criterion tries to approximate the set of separable 
states from the outside, our new criterion will work from 
the inside. Our criterion not only detects separability, 
but also provides an explicit separable decomposition of 
the separable states in terms of an integral over the Haar 
measure. Moreover, it has the same proven time com- 
plexity as the original DPS criterion. Due to its internal 
structure, the new criterion can be easily modified to 
perform linear optimizations over entanglement breaking 
channels. In particular, given any state estimation prob- 
lem, the method can output a sequence of measure-and- 
prepare strategies arbitrarily close to optimal. In the last 
pages of this letter, as an illustration of its efficiency, we 
will show its performance by computing the robustness 
of entanglement of a family of quantum states. 

But, first, some remarks about notation. In this let- 
ter, we will be mainly concerned with a finite dimensional 

def 

bipartite Hilbert space H = TLa^H-b- The set of all lin- 
ear operators acting on H will be denoted as B(H), and 
we will use the term state in order to refer to normal- 
ized non-negative operators in B{TL). Finally, the cone 
of separable operators, i.e., the conical combination of all 
pure product states a)^ a\ ® I^bX^bI}, will be called 
S. Note that any operator in S must be necessarily non- 
negative. 

The sequences of sets employed to approximate iS in 
the DPS criterion are either the sets {S^} of iV-(Bose) 
symmetrically extendible operators or the sets {5^} 
of PPT 7V-(Bose) symmetrically extendible operators. 
These sets are defined as follows: 

Aab € B(TLab) belongs to S N iff there exists an opera- 
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tor A ab n e B(Hab n ) that satisfies the following three 
conditions: 



f. A 



AB N 



> 0. 



2. tr B jv-i (A^jv) = Aab- 

3. A^n is Bose symmetric in H B N , i.e., A ab n(I a ® 
P^y m ) = A ab n, where Pjy m denotes the projector 



onto the symmetric subspace of 7i B 

Similarly, Aab G iff there exists an operator A^n S 
B(H ab n) fulfilling f-3 and the additional constraint: 



DPS 



lim^S K ,S»=S 



DPS* 




(lim_ v ^ .V • . v.; } .V 



4. A^n has a Positive Partial Transpose (i.e., it is 
PPT) [I w.r.t. thebipaxtition ABW 2 1|bW 2 J. 

Doherty et al. [3] proved that both sequences {iS^J^j 
and {Sp}^ =1 converge to S from the outside: 



FIG. 1: Two different ways of approximating the set <S of 
separable operators (dashed line): from the outside (DPS cri- 
terion) or from the inside (DPS* criterion). 



find a similar decomposition for any state pab & S„ m 



0- 



Note that, in both definitions, taking the limit N >> d 



The DPS criterion consists, precisely, in checking if implies, on one hand, that 5 w nS N ,S^ 
Pab £ SfL for all iV. From the relations above, the DPS 
criterion is clearly complete: if pab is entangled, then 
there exists an M such that pab Sfy , so all entangled 
states can be eventually detected. As all these sets are 
defined through linear matrix inequalities, the problem 
of determining whether a given state belongs to one of 
them can be cast as a semidefinite program (SDP) [l5j |. 

In [l6| it was shown that, for any pab G S n or Sp , 
a small perturbation on system B results in a separable 
state pab- More concretely, define d = dim(H B ) and 



ejv = 



{d-2,N mod 2) 



2{d - I) min 't 1 - T : P Lw/2J+i 



(.*) = 0}, (2) 



with Pn" 1 ' (x) being the Jacobi polynomials (20| . Then, 
the following closed cones of operators 



S N = { 



N 



N + d 



0~AB 



N + d 



a A ®l B :a AB eS"}, (3) 



S? = {(1 - e N )o AB + e N a A ® : a AB £ S£], (4) 

satisfy S N ,Sp C S, for all JV. Moreover, for any given 
Pab G we can find a corresponding separable de- 

composition in terms of an integral over a Haar measure. 
Suppose, for instance, that o~ ab n is a Bose symmetric 
extension of a AB G S n . Then, the separable decomposi- 
tion of the state pab G S n that results from perturbing 
(jab as in equation © is given by 



S^ [23], and, 

on the other hand, that Sp,S N w 5. It follows that 
the closures of the limiting sets limAr^oo S N ,Sp coincide 
with the set of separable states. 

Therefore, as opposed to the DPS criterion, the se- 
quences of sets {(S^Jjv, {Sp}iy converge to the set of sep- 
arable states from the inside. These sets can be charac- 
terized using semidefinite programming; in fact, a small 
modification of the computer codes employed to search 
over the sets S N , Sp allows to perform optimizations 
over the sets S N ,Sp. We thus arrive at a novel algo- 
rithm for entanglement detection complementary to the 
DPS criterion, and with a very similar structure 24|. In 
view to its strong link with the latter criterion, from now 
on we will refer to this new algorithm as DPS* . 

We will now proceed to study the speed of convergence 
of DPS*. First, notice that, as S C 5^,5^, equation (0J) 
tells us that we can perturb any general separable state 
Pab to a state p AB inside S N and Sp . This observation, 
together with the techniques employed in [lq . allows us 
to derive the following proposition: 



Proposition 1. Denote by | • || p the p-norm, i.e., 
\\Z\\ p = {tr(y/Z2^ P )) 1 /P, and let p £ S be a separable 
state. Then, there exists a normalized state a £ S N such 
that: 



2(d- 1) , „ d— 1 

\\p-<7\\i < AT | j , and ||p -alloc < 



N + d 



N + d' 



PAB = C J d4> B tr B N {l A (g) (f)% N a AB N } 



(5) 



where C is a normalization constant and d<j) B is a Haar 
measure over all pure states <p B on TLb- The reader can 



Analogously, there exists a normalized state cr p G Sp 
such that 

2(d~l)e N (d-l)e N 
\\p-a p \\i < , and \\p- < . 



N 
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It is worth noting that all the above bounds cannot be 
improved, since they correspond to the exact distances 
of Pab to the sets S N or whenever pab is a pure 
product state. 

Proposition [T] allows us to study the efficiency of the 
criterion based on the sequence of sets (S N ) as opposed 
to the one based on (S^ ) when applied to solve the weak 
membership problem of separability (WMEM(S)) 0,0, 
i.e., the problem of determining whether a state pab is 
separable or not up to some precision S. Since our bounds 
on the speed of convergence of DPS* have the same scal- 
ing as those derived for the DPS criterion in [lij ]. our 
conclusions on the complexity of the former cannot but 
be the same, namely: 

1. From the point of view of time complexity, it is 
always preferable to consider the sets rather 
than S . Whereas the time complexity of the lat- 
ter scales with S as O (d A (ki/5) 6dB ) , the number 
of operations required under the PPT constraint 
scales as O (d\{k 2 /5) 4dB ). 

2. From the point of view of space complexity, some- 
times it may be more convenient to search over the 
sets {i^} rather than over {5^}. 

3. In any case, due to its fast proven convergence as 
compared to other methods (|, the DPS* criterion 
is one of the most efficient algorithms for entangle- 
ment detection. Because of the polynomial depen- 
dence on <1a of its bounds on time or space com- 
plexity, DPS* is specially useful to attack problems 
where the dimension of one of the systems is much 
bigger than the other one's. 

We will now illustrate the power and versatility of this 
algorithm by showing how to use it to obtain approxi- 
mate solutions for two important problems that appear 
frequently in quantum information: 
a) State Estimation Problems 

In a generic pure state estimation scenario, a source 
randomly chooses a pure state out of a probability 
distribution pi and then encodes it into a state ^ to 
which we are given full access. The state estimation prob- 
lem consists on finding the optimal measure-and-prepare 
strategy that allows us to reconstruct the original state 
with high fidelity. In [IS], it is explained how to map 
state estimation (SE) problem problem into a linear op- 
timization over the set S of separable operators, via the 
relation 



F = max{tr(pAsAAs) : A A b e S, A A = I}.. 



(0) 



where F is the optimal average fidelity and pab = 
YliEi^'i ® is given by the particular SE problem. 
In [19( it is also shown that any separable decomposi- 
tion Aab = Y^x M x <8> <j) x of an operator Aab satisfying 



the above conditions can be interpreted as a measure- 
and-prepare strategy consisting of applying the POVM 
{Mx}x and preparing the state <j> x depending on the out- 
come x. 

Now, consider the sequence of optimization problems: 



■ V _ 



max{tr(p^ B A4 S ) : A A b € S , A A = I}, 
m&x{tii P ab A AB ) ■ A AB € A A = I}. (7) 



These optimizations can also be cast as a semidefinite 
program. Therefore, they can be efficiently computed as 
long as the index N is not very high (note nevertheless, 
that, in general, if we fix N for increasing ds, we will 
end up with a very bad approximation). Also, because 
we are optimizing over particular regions of S, in general 
we will get a suboptimal result, i.e., F N , F^ < F for all 
N. However, it is clear that limjv^oo F N , F^ = F. 

As we have already mentioned, we can write down a 
separable decomposition for the operator Aab output by 
the computer, via E q.Q in case Aab € S N or by means 
of a similar formula [16| in case Aab S Sp • This separa- 
ble decomposition can be subsequently interpreted as a 
measure-and-prepare strategy; for any SE problem, the 
DPS* method can thus provide us with a sequence of 
state estimation strategies that converge asymptotically 
to the optimal one. 

It would be natural to wonder how fast this conver- 
gence is, i.e., how far F N is from F for finite N. Follow- 
ing the lines of (l6j , we arrive at a sequence of upper and 
lower bounds on F given by 



N 



F N < F < F 



F N <F<F N 
p — — p 



N 



1 



T^( P P N 
1 - e N » 



(8) 



To get an idea of the efficiency of these algorithms, we 
refer the reader to [3, since the lower bounds on the 
maximal fidelity of the state estimation problems con- 
sidered in that paper actually correspond to the values 
computed through DPS*. 

b ) Computation of the robustness of entanglement 

Let p £ B(H) be a quantum state. We will define its 
(separable) robustness of entanglement [13] R(p) as 



R(p) = min{tr(er) : er G <S, p + cG S}. 



(9) 



Used in combination with the DPS criterion, the DPS* 
criterion allows us to determine the robustness of en- 
tanglement of any quantum state pab up to arbitrary 
precision. Indeed, note that, optimizing over S N (Sp) 

or S N (Sp) instead of iS in the previous definition, wc 
would obtain lower and upper bounds on R(p), respec- 
tively. And, of course, both optimizations can be per- 
formed using semidefinite programming. 

In [l8| , the authors introduced a class of uniparametric 
families of n-partite quantum states. Given a unitary 



4 



operator V acting over n — 1 qubits, the corresponding 
family of states is denned as 



i 



i 



i 



Ii ® I23...T. + ™ a|0)(l| ® V + — a|l)(0| ® U* 



The separability properties of these states are very im- 
portant, for if one could prove that any such state p^ 
is multiseparable for all a < l/poly(n), then one would 
have an example of a quantum computation that super- 
sedes any classical algorithm but nevertheless does not re- 
quire entanglement [181 ]. For the case n = 3, the authors 
showed that, for any unitary V, the state p% is always 
separable with respect to the partition 1|23. Moreover, 
their numerical tests suggested that p\ is PPT with re- 
spect to the partition 1 2 1 3 for all a < 1/2. Nevertheless, 
they conjectured that, in some cases, the state should 
remain entangled for lower values of a. The DPS* cri- 
terion presented above strongly suggests that this is not 
the case. 

We generated 1000 random unitaries {Vi]i C SU(A) 
according to the Haar measure and applied the DPS* cri- 
terion to derive upper bounds on the robustness of entan- 
glement of the corresponding states {pX%}i by consider- 
ing extensions over the last qubit. We used the MATLAB 
package YALMIP [2l| in combination with SeDuMi [22 1 
to perform the numerical calculations. After an appro- 
priate optimization over the set 5 3 , all the corresponding 
upper bounds turned out to be zero, thus proving the 
separability of the previous sample of states. 

In order to discard statistical effects, we considered the 
unitary operator U — 2P|y m — 12 <8> I2, whose associated 
states p^ appear to have the grea test negativity of the 
whole family for fixed a > 0.5 |25j. Fig. [2] shows upper 
(DPS*) and lower (DPS) bounds on the robustness of 
entanglement for different values of a. This time, we 
optimized over the sets 6> 3 , S 15 and Sp, respectively. It 
is clear that for all values of a below 0.5, the state p„ is 
separable. 

In conclusion, in this letter we have introduced a new 
criterion for separability detection. This criterion has 
been inspired by and it is in a sense complementary to the 
one conceived by Doherty et al. [l|. Whereas the latter 
one aims at approximating the set of separable states 
from the outside, our new method works from the inside, 
i.e., by defining a sequence of sets of states contained in 
S. From pj}, it follows that the method can be easily 
extended to deal with multiseparability problems. 

This new criterion works basically by taking the states 
defined by Doherty et al. and applying a perturbation 
to make them separable. We believe, however, that the 
size of this perturbation is larger than required: in the 
PPT case, at least for small dimensions of Alice's sys- 
tem, a smaller transformation (in some cases the iden- 
tity) should be enough to guarantee the separability of 
the output state. Note that the bounds on the distance 
to arbitrary separable states are actually independent of 




FIG. 2: Upper (DPS*) and lower (DPS) bounds on the en- 
tanglement robustness R of as a function of a. The op- 
timizations over the sets <S 3 , <S 15 and Sp correspond to the 
dashed-dotted, dashed and solid lines, respectively. 



d,A- If future research found an optimal linear map to 
turn Doherty et al.'s states into separable states that 
took into account the dimensionalities of both systems, 
the speed of convergence of the resulting improved DPS* 
criterion would be much faster. 

Nevertheless, a wide range of applications follow from 
our method in its present state. We have seen that the 
method alone can be used to determine the best experi- 
mental setup in quantum tomography protocols, and we 
have also shown how the combination of both the DPS 
and DPS* methods allows to solve with arbitrary pre- 
cision computationally hard problems, like the calcula- 
tion of the entanglement robustness. Since some impor- 
tant problems in Complexity Theory like PARTITION 
or CLIQUE can be reduced to the separability problem 
it is not unrealistic to expect that, with time, our 
method will find new applications outside the scope of 
Quantum Information Science. 
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EPSRC QIP-IRC and is supported by EPSRC grant EP 
C546237 1, the Royal Society, the EU Integrated Project 
QAP and an Institute for Mathematical Sciences postdoc 
fellowship. 
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for N » 1 \2C 



where 



(d-l)N' A 

jd-2,1 denotes the first positive zero of the Bessel func- 
tion Jd-2(x). It follows that limjv^oo cm = 0. 
In this case, though, we do not have a hierarchy anymore, 
i.e., in general, S N <£ S N+1 . 
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